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Abstract 

A  simple  model  based  on  the  mean  spherical  approximation  is  described  for  the  dielectric 
properties  of  polar  solvents.  The  molecular  properties  required  arc  the  dipole  moment,  the 
molecular  polarizability,  and  the  molecular  diameter.  It  is  shown  that  inclusion  of  a  stickiness 
parameter  accounting  for  other  electrostatic  features  and  non- sphericity  in  the  real  system  leads  to  a 
simple  analytical  result  from  which  the  dielectric  constant  may  be  calculated.  Th*stickiness 
parameter  is  shown  to  be  correlated  to  the  magnitude  of  dipole-dipole  interactions  in  the  system. 


Introduction 

Application  of  the  mean  spherical  approximation  (MSA)  to  the  description  of  the  dielectric 
properties  of  polar  solvents  dates  from  the  work  of  Wertheim.^  In  this  model,  solvent  molecules 
are  represented  as  hard  spheres  with  centrally  located  point  dipoles.  When  applied  to  water,  the 
estimated  dielectric  constant  is  considerably  less  than  the  experimental  value.^  Thus,  one  must 
replace  the  experimental  dipole  moment  by  a  higher  effective  value  in  order  to  obtain  the 
experimental  dielectric  constant  from  the  MSA  model.^*^  It  is  generally  recognized  that  the  failure 
of  models  for  water  based  on  a  hard  sphere  -  point  dipole  description  is  due  to  the  fact  that  this 
simple  molecule  has  a  more  complex  multipole  character.**  The  complex  electrostatic  description 
leads  to  its  well  known  hydrogen-bonded  character. 

Application  of  the  MSA  model  to  larger  polar  solvents  is  certainly  expected  to  result  in  other 
problems.  Thus,  if  one  considers  the  series  of  normal  alcohols,  the  representation  of  the  molecule 
as  a  sphere  with  a  central  point  dipole  becomes  worse  as  the  length  of  the  alkyl  chain  increases. 
Problems  of  a  similar  nature  are  encountered  when  one  considers  a  series  of  aprotic  solvents  such 
as  acetonitrile,  and  alkyl  cyanides  with  larger  alkyl  groups.  Thus,  in  general,  one  recognizes  two 
problems  with  any  model  for  a  polar  liquid  based  on  hard  spheres  with  centrally  located  point 
dipoles.  The  frrst  is  that  the  electrostatic  description  of  the  molecule  may  require  that  the 
quadrupole  moment  or  higher  order  multipoles  be  considered.  The  second  problem  is  that  the 
molecule  may  not  be  well  represented  as  a  sphere. 

One  way  of  dealing  with  the  above  inadequacies  is  to  introduce  stickiness  in  the  model.  In 
this  way,  one  can  correct  for  the  multipole  nature  of  the  polar  molecule  and  its  soft  intermolecular 
interactions.  The  sticky  interaction  either  aligns  or  misaligns  the  electrical  dipoles  of  the  solvent 
molecules,  and  therefore  can  be  called  a  structure  forming  or  structure  breaking  parameter.  In  the 
present  paper,  the  original  work  of  Wertheim'  is  extended  by  including  a  sticky  interaction 
parameter  in  the  description  of  the  polar  molecule.  Otherwise,  the  solvent  molecule  is  considered 
to  be  a  hard  sphere  possessing  a  permanent  point  dipole  and  a  polarizability.  The  model  is  applied 
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to  representative  polar  solvents,  including  both  aprotic  and  protic  systems,  in  order  to  assess  the 
role  of  the  stickiness  parameter  in  determining  bulk  dielectric  properties. 

The  Model 

On  the  basis  of  Wenheim's  original  work^,  one  defines  a  polarization  parameter,  X,  by  the 
relationship 

X2(X+l)4=16es  (1) 

where  £5  is  the  static  dielectric  constant  of  the  polar  liquid.  The  parameter,  X,  also  arises  when  the 
MSA  is  applied  to  the  estimation  of  ion-dipole  interactions  in  an  infinitely  dilute  electrolyte 
solution^,  and  is  designed  to  take  into  account  softness  and  nonspherecity  of  the  solvent  core.  X  is 
obtained  from  the  energy  parameter,  b2,  which  comes  from  the  solution  of  the  MSA  for  dipole- 
dipole  interactions  for  a  system  composed  of  hard  spheres  with  centrally  located  point  dipoles. 

The  relationship  between  X  and  b2  is 


1 

^  =  l-b2/6 


(2) 


where 

00 

^2  =  J  dr^-  ---  Gs^  (3) 

o 

A 

h'  '2  (r)  is  the  relevant  coefficient  from  the  invariant  expansion  of  the  total  pair  correlation  function 
h(r)  which  depends  on  the  distance,  r,  between  two  dipolar  spheres.^  Ps  is  the  number  density  of 
molecules  in  the  liquid  and  Os.  their  diameter.  The  excess  internal  energy  per  unit  volume  for  the 
liquid  is  given  by 


£  d2  b2 


(4) 


where 


2  47t  Ps  p2 

2  ■  3kBT 


(5) 
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Here,  p,  is  the  permanent  dipole  moment  of  the  solvent  molecule,  and  ke,  the  Boltzmann  constant. 
In  the  MSA,  when  the  molecules  arc  non-polarizable,  one  obtains  the  result 


,2  (X+2)2 

= - 9 - 


■-t) 


(6) 


Wenheim  also  studied  the  MSA  and  other  approximations  for  the  case  of  dipole  hard  spheres 
which  arc  polarizable.^  This  work  was  reformulated  to  show  that  it  is  equivalent  to  a  mixture  of 
hard  spheres  of  equal  diameter  but  varying  dipole  moment.^  As  a  result,  when  the  molecules  have 
a  polarizability  a,  eq.  (6)  becomes 


o  (X+2)2 
= - 9 - 


1-1 

Es 


(7) 


where 


2  2 

dp  =  a2  d2  +  4n  a  ps  a 


with 


a-i  =  1  - 


4  b2  g 
3os^ 


(8) 


(9) 


Solving  eq.  (7)  for  the  parameter  X,  one  obtains 

X  =  (1 +3da)i/2  .  1  (10) 

where 

da  =  dp/(l-l/£s)i/2  (11) 

When  one  applies  the  above  system  of  equations  to  water  for  which  p  is  equal  to  1.83 
Debyes  and  ai,  to  1.47x10’^  nm^,  one  obtains  the  result  that  dp  =  4.33  and  X  =  2.75.  As  a  result, 
on  the  basis  of  eq.  (1),  the  MSA  estimate  of  Es  for  water  at  25  °C  is  93.7,  which  can  be  compared 
with  the  experimental  result  of  78.3.  This  suggests  that  a  small  correction  to  the  contact  correlation 
functions  would  yield  good  agreement  with  experiment.  This  is  achieved  here  using  the  directional 


sticky  model  described  by  Blum,  Cummings  and  Bratko.^  In  its  simplest  dipolar  version, 
introducing  a  longitudinal  stickiness  parameter  to,  eq.  (7)  becomes 
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Having  estimated 


_,2  x2(X+2)2  r..  .  n 

= - 9 -  [d-  toA)2  -  ^  (12) 

on  the  basis  of  eq.  (8),  one  can  then  determine  the  stickiness  parameter  to  by 


solving  eq.  (12)  using  the  value  of  X  obtained  from  the  experimental  value  of  £5  (eq.  (1)).  One 
then  has  a  self-consistent  model  in  which  the  experimental  value  of  e®  is  estimated  from  the 

molecular  dipole  moment  and  polarizability  using  the  experimentally  defined  polarization  parameter 
X  with  an  additional  stickiness  parameter  to. 


This  model  is  now  examined  for  selected  polar  solvents.  However,  before  examining  the 
static  solvent  permittivity,  the  high  frequency  permittivity,  and  the  ability  of  the  MSA  to  estimate 
the  refractive  index  from  the  molecular  polarizability  are  considered. 

Solvent  Permittivity  at  Optical  Frequencies 

At  optical  frequencies,  the  system  only  responds  to  an  electrical  field  via  its  molecular 
polarizability.  Thus,  in  this  case  the  parameter  dp  is  given  by  the  equation 

dp  =  ^napstt  (13) 

Furthermore,  the  polarization  parameter,  X,  is  now  calculated  from  the  optical  permittivity,  Eop, 
that  is,  from  the  equation 


Eop  =  >.2(U1)4/16  (14) 

where  Eop  is  the  square  of  the  refractive  index.  Thus,  an  alternate  way  to  estimating  dp  is  via  the 
equation 


2  _  X2(X-^2)2 
P  *  9 


(15) 


The  parameters  b2  and  a  defined  above  by  equations  (2)  and  (9),  respectively,  are  now  estimated 
using  the  value  of  X  obtained  from  eq.  (14). 


In  the  case  of  water  at  25  °C,  the  valuer  of  Eop  is  1.7756.  The  corresponding  high  frequency 
value  of  the  polarization  parameter  X.  is  1.151  on  the  basis  of  eq.  (14).  One  may  now  estimate  the 
parameter  dp  by  using  eq.  (15)  and  then  calculate  the  molecular  polarizability  on  the  basis  of  eq. 
(13).  For  water,  the  MSA  estimate  of  a  is  1.483  x  10'^  nm^,  a  result  which  compares  very  well 
with  the  value,  1.47 1  x  lO^  nm^,  estimated  on  the  basis  of  the  Loren tz- Lorenz  relationship.^  In 
the  latter  case,  the  polarizability  is  calculated  from  the  molar  refraction  R  defined  by  the  equation 

eop+2 

w  here  ^  is  the  molar  volume  and  Nq,  the  Avogadro  constant.  Keeping  in  mind  the  fact  that  the 
value  of  the  polarizability  depends  very  much  on  whether  it  is  obtained  from  gas  or  liquid  phase 
measurements it  seems  appropriate  that  it  be  obtained  from  the  refractive  index  of  the  liquid. 

Values  of  a  were  calculated  for  all  other  solvents  considered  using  both  the  MSA  and  the 
Lorentz-Lorenz  relationship.  In  all  cases,  the  MSA  estimate  is  slightly  higher,  typically,  by 
approximately  1%.  Because  the  difference  is  so  small,  the  MSA  result  was  used  in  all  further 
calculations  and  is  tabulated  in  Table  1  together  with  other  pertinent  solvent  data. 

The  Static  Solvent  Permittivity 

The  present  model  was  applied  to  20  different  polar  solvents,  both  protic  and  aptotic,  chosen 
so  that  the  role  of  the  stickiness  parameter  could  be  assessed.  Three  solvent  groups  considered 
w  ere  water  and  the  lower  alcohols  in  which  the  polar  group  is  -OH,  the  alkyl  cyanides  with  -C  s  N 
as  the  polar  group,  and  the  amides  whose  dielectric  parameters  depend  markedly  on  whether  the 
solvent  is  ptodc  or  aprotic.  The  dielectric  and  other  relevant  parameters  for  these  systems  are 
summarized  in  Table  1.  Some  comment  should  be  made  about  the  values  of  the  dipole  nwment 
chosen  from  the  compilations  of  McClellan.  *2  The  value  of  this  parameter  depends  somewhat  on 
whether  it  was  measured  in  the  gas  phase  or  in  a  dilute  solution  in  a  non-polar  solvent.  The  chosen 
values  for  a  group  of  solvents  were  obtained  by  one  method,  preferably  using  measurements  from 


one  '.aboratory.  The  results  obtained  in  the  present  study  depend  markedly  on  the  value  chosen  for 
this  parameter,  so  it  is  important  that  it  properly  reflect  trends  within  a  given  group  of  solvents. 

The  stickiness  parameter,  to,  also  recorded  in  Table  1  was  estimated  from  the  dielectric  data 
in  the  following  way.  Initially,  the  polarization  parameter,  X,  and  the  related  energy  parameter,  b2, 
were  calculated  from  the  bulk  dielectric  constant,  £5,  (eqs.  (1)  and  (2)).  Then  the  molecular 
parameter,  di,  was  estimated  on  the  basis  of  eqs.  (5),  (8)  and  (9)  using  the  molecular  polarizability 
found  earlier,  and  the  dipole  moment,  density,  and  molecular  diameter  tabulated  in  Table  1. 

Finally,  the  value  of  to  was  obtained  by  solving  eq.  (12)  using  the  above  values  of  X  and  dp. 

From  the  results,  it  is  apparent  that  to  varies  over  a  wide  range  from  a  low  of  -1.498  for 
formamide  to  a  high  of  0.495  for  1 -propanol.  Negative  values  of  to  indicate  that  dp  is 
underestimated  when  it  is  calculated  on  the  basis  of  eq.  (8).  It  also  follows  that  the  MSA  then 
underestimates  the  ordering  effect  of  a  given  dipole  on  its  neighbours.  In  the  case  of  protic 
solvents  such  as  water  and  formamide,  this  is  attributed  to  the  presence  of  hydrogen  bonding,  and 
important  multipole  contributions  to  the  electrostatic  description  of  the  molecule.  Negative  values 
of  to  are  also  found  for  all  aprotic  solvents  considered.  This  also  reflects  the  fact  that  the  dipole 
moment  cannot  be  considered  to  be  centrally  located  in  any  of  the  polar  solvent  molecules 
considered. 

It  is  interesting  to  follow  the  change  in  to  with  the  strength  of  the  dipole-dipole  interactions  in 
the  system.  This  can  be  seen  for  the  water-alcohol  series  and  also  the  alkyl  nitrile  series.  Plots  of 
to  against  p^/d^  are  linear  for  both  series  over  a  range  of  4  kT  (Fig.  1).  What  is  interesting  about 
these  results  is  that  to  becomes  more  positive  as  dipole-dipole  interactions  decrease  in  magnitude. 
This,  in  turn,  can  be  related  to  the  dipolar  group  in  the  solvent  molecule  moving  further  and  further 
away  from  the  molecular  center.  No  such  trend  can  be  seen  for  the  amides.  Undoubtedly,  this  is 
due  to  the  fact  that  the  group  of  solvents  considered  includes  both  protic  and  aprotic  systems  so 
that  the  complete  electrostatic  description  varies  considerably  among  the  members  of  this  group. 
N-methylformamide  is  the  only  other  solvent  listed  which  has  a  positive  value  of  to.  Both 
formamide  and  N-methylformamide  are  systems  which  arc  extensively  hydrogen  bonded  but  with 
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quite  different  structures  as  a  result  of  change  in  number  of  acidic  hydrogens.  This  change 
undoubtedly  leads  to  the  large  difference  in  to  for  these  solvents. 

The  remaining  solvents  listed  in  Table  I  are  all  aprotic  and  possess  negative  values  of  to. 

This  parameter  is  plotted  against  the  dipole-dipole  interaction  energy  unit,  p^/d^  in  Fig.  2.  There  is 
clearly  a  correlation  between  these  two  quantities,  the  value  of  to  becoming  more  negative  as  p^/d^ 
increases.  This  suggests  that  the  smallest  molecules  with  largest  dipole  moments  depart  most 
severely  from  the  simple  MSA  result.  This  result  reflects  the  fact  that  the  tendency  for  local  order 
is  underestimated  by  the  MSA.  The  scatter  is  much  larger  than  that  seen  in  Fig.  1  where  molecules 
with  the  same  polar  group  were  considered.  Not  only  are  the  polar  groups  different  but  the  size 
and  shape  of  the  non-polar  pan  of  the  molecule  varies  considerably  among  the  solvents  considered. 
In  addition,  systems  like  nitrobenzene  and  benzonicrile  have  very  polarizable  phenyl  rings  whereas 
the  other  molecules  do  not.  In  conclusion,  the  observed  correlation  is  chiefly  a  result  of  the  change 
in  molecular  size,  and  corrects  the  MSA  model  for  the  fact  that  the  dipolar  group  is  not  centrally 
located  in  the  solvent  molecule. 

Discussion 

The  MSA  model  for  the  dielectric  propenies  of  polar  solvents  described  here  is  remarkably 
simple.  Using  the  molecular  dipole  moment,  polarizability,  and  diameter  estimated  from  gas 
solubilities,  one  is  able  to  calculate  the  relative  permittivity  with  two  additional  parameters,  namely, 
the  polarization  parameter,  A,,  and  the  stickiness  parameter,  tg.  The  first  of  these,  X,  is  calculated 
directly  from  the  static  dielectric  constant.  Furthermore,  the  stickiness  parameter,  t©,  is  clearly 
correlated  to  the  ratio  which  is  proportional  to  the  dipole-dipole  interaction  energy.  As  the 
size  of  the  non-polar  group  increases  in  a  series  of  molecules  with  the  same  polar  group,  the  ratio 
pVd?  decreases,  and  tg,  increases.  The  correlation  between  these  quantities  is  perfectly  linear  for  a 
given  series  of  compounds. 

The  limitations  of  the  simple  MSA  in  describing  the  properties  of  polar  liquids  are  well 
known  the  present  extension  adding  very  little  to  the  complexity  of  the  model.  If  one  compares 
the  results  of  the  present  model  to  those  from  molecular  dynamics  calculations,  it  is  clear  that  the 
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latter  approach  only  would  yield  reasonable  results  if  a  detailed  electrostatic  description  of  the 
solvent  molecule  is  provided.  The  fact  that  the  MSA  provides  an  analytical  model  with  empirical 
parameters  which  can  be  qualitatively  related  to  known  properties  of  the  system  suggests  that  it  is  a 
very  good  basis  for  describing  the  dielectric  properties  of  polar  liquids. 
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Legends  for  Figures 


Figure  1. 


Figure  2. 


Plot  of  the  stickiness  parameter,  to,  against  the  dipole-dipole  interaction  energy  factor 
p2/d3  for  the  lower  alcohols  and  water,  and  for  the  alkyl  cyanides.  The  abbreviations 
for  the  solvent  given  next  to  each  {XDint  are  defined  in  Table  1. 

As  in  Figure  I,  but  for  all  aprotic  solvents  listed  in  Table  1. 
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